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Abstract 

The problem of the origin of the /i parameter in the Minimal Supersymmetric Standard 
Model can be solved by introducing singlet supermultiplets with non-renormalizable 
couplings to the ordinary Higgs supermultiplets. The Peccei-Quinn symmetry is bro- 
ken at a scale which is the geometric mean between the weak scale and the Planck 
scale, yielding a u term of the right order of magnitude and an invisible axion. These 
models also predict one or more singlet fermions which have electroweak-scale masses 
and suppressed couplings to MSSM states. I consider the case that such a singlet 
fermion, containing the axino as an admixture, is the lightest supersymmetric parti- 
cle. I work out the relevant couplings in several of the simplest models of this type, 
and compute the partial decay widths of the next-to-lightest supersymmetric particle 
involving leptons or jets. Although these decays will have an average proper decay 
length which is most likely much larger than a typical collider detector, they can occa- 
sionally occur within the detector, providing a striking signal. With a large sample of 
supersymmetric events, there will be an opportunity to observe these decays, and so 
gain direct information about physics at very high energy scales. 



Contents 

1 Introduction 



2 Singling masses and couplings in models with an intermediate-scale solu- 



tion to the \i problemj 7 



3 Mixing of the singling with MSSM neutralinos| 10 



4 Representative results for decays to the singling] 13 



4.1 Neutralino decays! 13 



4.2 Slepton decays] 17 



5 Conclusions! 18 



Appendix: Complete formulas for neutralino decay widths including the 



effects of arbitrary phases] 19 



1 Introduction 

Supersymmetry (for reviews, see |l|, §]) has the ability to stabilize the hierarchy between 
the electroweak and Planck scales. However, the minimal supersymmetric standard model 
(MSSM) still requires an explanation for the magnitude of the supersymmetric Higgs mass 
parameter /i. Assuming that there are no fine-tuned cancellations in the MSSM Higgs 
potential, /i should be of roughly the same magnitude as the soft supersymmetry-breaking 
masses. This suggests that /i arises as a vacuum expectation value (VEV) which is fixed 
by a potential with dimensionful parameters that are in turn determined by supersymmetry 
breaking. 

Supersymmetry also requires some additional structure in order to solve the strong CP- 
problem. The fi parameter breaks the Peccei-Quinn (PQ) symmetry |J that is otherwise 
naturally present in the MSSM at the renormalizable level, so it is an attractive proposition 
that the dynamics which leads to the /x term simultaneously provide for an invisible axion. 
Astrophysical constraints on the axion leave open a window ||] from roughly 

10 9 GeV < / < 10 12 GeV (1.1) 

for the VEV of the PQ-breaking field. 
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In this paper I will consider the phenomenology of a class of models with an invisible axion 
P, |6j of the DFSZ type || [7], |8j. In these models, the \x term arises from non-renormalizable 
terms in the superpotential, for example: 



(Here Mp = 2.4 x 10 18 GeV is the reduced Planck mass, and is a dimensionless coupling 
which I assume is not much smaller than unity.) The sum of the PQ charges of X± and X 2 
must be equal and opposite to that of the MSSM Higgs superfields. When the scalar compo- 
nents of the neutral chiral superfields Xi and X 2 acquire VEVs of order /, the approximate 
PQ symmetry is spontaneously broken, giving rise to a fi term of the right order of magni- 
tude and an invisible axion. It is natural to assume that this occurs with X\ and X 2 along a 
nearly flat direction in the potential. Then the low-energy degrees of freedom will typically 
include a pair of neutral chiral supermultiplets which are mixtures of the original Xi and X 2 . 
One of these contains as its imaginary scalar component the invisible axion of the model, 
and its fermion superpartner, the "axino" . In some models, X\ and X 2 are the same field, 
so that the axino is the only new light singlet fermion. If X\ and X 2 are distinct, then there 
will be another light Majorana fermion "singlino". The axino and singlino both have odd 
i?-parity and can mix. They obtain masses which are not much larger than the weak scale 
(but might be as small as of order a keV depending on the details of the model ||). The 
upper limit can be understood from the facts that the (mass) 2 splittings between members 
of the same supermultiplet are bounded above by roughly m 2 .^, the squared gravitino mass, 
and the axion is nearly massless. 

The axino and the singlino are very weakly coupled to MSSM states, and cannot be 
directly produced in collider experiments at any significant rate. However, if either (or both) 
of these particles is lighter than all of the MSSM superpartners, then it will be the light- 
est supersymmetric particle (LSP) and can appear in decays from ordinary supersymmetric 
events.[] In this paper I will argue there will be an opportunity in the era after supersym- 
metry is discovered to search for and measure the very long lifetime of the next-to-lightest 
supersymmetric particle (NLSP) into final states that include the singlino or axino, despite 
its very weak coupling. 

In order to discuss the phenomenology in a general way, I will denote the relevant axino 

'If i?-parity is conserved and the singlino or the axino is the LSP, it will be absolutely stable and 
could dominate the energy density of the universe too soon. This potential problem can be solved (as in 
many similar cases) by invoking a low reheat temperature, at the cost of requiring in addition a low-scale 
baryogenesis mechanism. Furthermore, the NLSP decays will safely occur long before nucleosynthesis in the 
standard cosmology. 



w 



Mp 



XiX 2 H u Hd- 
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or singlino LSP by S, and refer to it generically as a singlino. It is part of a superfield S. In 
the low-energy theory it participates in a superpotential of the form 



W = fi\l + -Sj H u H d + -m^S 2 (1.3) 

in which 

H = -^(X 1 )(X 3 ), (1.4) 

and I have introduced a dimensionless coupling parameter 

e~v/f (1.5) 

with v = 175 GeV, the electroweak scale. (There are also soft mass terms for the scalar 
components of S, which will not concern us.) For example, if X\ and X 2 are the same field, 
then one can read off from eq. ( |1.2| ) that e — 2v/f. For numerical purposes, this paper will 
use as a benchmark the value e m 10~ 8 corresponding to / ~ (few)xlO 10 GeV. There is also 
a dimensionless, holomorphic soft term in the lagrangian 

Hh 

— ^SUSY breaking = ~7~7 X\ X 2 H u Hd + C.C. (1.6) 

1V1 p 

where hj, is of order mw- This gives rise to (among other terms) the necessary holomorphic 
soft (mass) 2 term for the Higgs bosons in the MSSM: 

— ^SUSY breaking = bH u Hd + C.C. (1.7) 

Note that b = h b (X 1 )(X 2 )/M P is of order ?77 ^/ j clS required for proper electroweak symmetry 
breaking. 

The coupling e and mass parameter parameterize our ignorance of the high-energy 
theory. The smallness of e means that S nearly decouples. However, the conservation of 
.R-parity implies that if the singlino is the LSP, then decays of the NLSP to S will not 
suffer any competition and can be observed if they happen within a collider detector. These 
decays occur and are potentially observable because the "singlino" S mixes slightly with 
the gauginos and Higgsinos, as well as couples directly to higgsino-Higgs pairs. In that 
sense, these models are similar to the well-studied [[ITJ-[[nj "next-to-minimal supersymmetric 
standard model" (NMSSM) ||16|| . The differences include: the extremely small magnitude 
of e; the fact that the field S does not obtain a VEV; the absence (or at least weak-scale 
phenomenological irrelevance) of an S 3 term in the superpotential; and the presence of a tree- 
level supersymmetric mass term for S. Nevertheless, it is useful to compare the situation 
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under study here to a very weakly coupled limit of the NMSSM. Indeed, the possibility of 
macroscopic decays involving a singlino have already been noted in ref. [pi, but considering 
larger couplings (effectively e > 10~ 6 ) and for smaller values of rrig appropriate for LEP. 

If the NLSP is the lightest of the ordinary MSSM neutralinos Ni, then it can decay 
according to 

Ni - fJS (1.8) 

through virtual sleptons and squarks and virtual or on-shell Z bosons and Higgs bosons. 
The decay width is estimated very roughly by 

m N 

T~ L |e/i/f| x (suppression factors). (1.9) 

167T 

The suppression factors include the effects of electroweak couplings, mixing angles, kinematic 
suppressions, and (if the mediating boson is not on-shell) three-body phase space. Without 
these effects, the rough estimate (for ~ 100 GeV and |e/i/t>| ~ 10 -8 ) would be of order 
1 meter -1 for T. 

After including the suppression effects in realistic models, we will find that when N\ is 
allowed to decay through an on-shell CP-even Higgs boson h°, the inverse decay width is of 
order meters or tens of meters. Of course, the decay Aq — > h°S may not be kinematically 
allowed. In that case, there may still be allowed decays Aq — ► Z°S. These are typically fur- 
ther suppressed by a mixing angle, because the singlino must mix with the MSSM higgsinos 
in order to couple to the Z boson. Nevertheless, we will find that these mixing angles are 
typically large enough so that the inverse decay widths can be of order hundreds of meters. 
Finally, it may be that — < mz- In that case, there can still be three-body decays 
Aq — > ffS through virtual sleptons, squarks, and the Z°. (Decays through off-shell Higgs 
bosons can also occur, but are typically very small because the MSSM Higgs boson widths 
are tiny unless they are heavy.) If the Z boson is far off-shell, then with the usual model 
prejudices that sleptons are much lighter than squarks one finds that the smallest inverse 
decay widths are for £ + £~S final states, and can be of order tens of kilometers. (All of these 
results assume e ~ 10~ 8 , and the decay widths must be scaled with e 2 .) 

The majority of decays Aq — > ffS will evidently occur well outside of a typical collider 
detector. However, with a significant number of supersymmetric events available, a small 
but finite fraction will occur inside the detector where the displaced secondary vertex can be 
distinguished. Since the decaying Aq and the resulting S are invisible, the experimental sig- 
nature will involve an energetic lepton-antilepton pair or dijet pair with a significant opening 
angle appearing "from nothing" (with no corresponding charged particle track pointing back 
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to the interaction point) at a common point. This determination could be accomplished 
within the inner tracking volume of a detector, but might also be possible and perhaps even 
easier to distinguish if the decay occurs farther from the beam pipe. Thus, decays occurring 
within a meter to several meters from the interaction point could be a striking, if rare, signal. 

At the Large Hadron Collider, the number of supersymmetric events expected per year 
with a low luminosity option of 10 fb _1 /year can be roughly of order a few thousand to 
a few million or more for 200 GeV > > 50 GeV. (This assumes m~ ~ m~ w 7m^ ; 
of course this is quite model-dependent.) Every event gives two possible NLSP decays. 
Therefore one can aspire to detect rare decays with widths as long as hundreds of kilometers 
by searching within the supersymmetric event sample fl7fl . In the limit of small decay 



widths, the probability that a particular N\ — > ffS decay will occur within a distance L of 
the interaction point is given by 

P(L)^LT/P 7 (1.10) 

where T is the invariant partial width for that decay channel. Supersymmetric events will 
be "tagged" by the other particles from the sparticle decays and the presence of large $ T . 
Since the slowly decaying N\ will be massive and not ultrarelativistic, one can use timing 
information together with the pointing information from the tracking detectors and drift 
chambers and perhaps vetoes from the muon system to eradicate backgrounds from cosmic 
rays and other sources. There have also been proposals motivated by gauge-mediated su- 
persymmetry breaking models and by neutralino decays to axino and photon [18] to build 



special detector components and instrumented tunnels to aid in the search for very slow 
decays JT^, |2"0|| . At future e + e~ linear colliders, supersymmetric event rates are smaller, but 
the two-body decays mentioned above might occur often enough to be detected. The signal 
is somewhat more problematic at future runs of the Fermilab Tevatron collider, since the 
total supersymmetric event rates are likely to be considerably smaller. In this paper, I will 
simply remain optimistic and choose to present results for decay partial widths down to as 
small as (1000 km)" 1 . 

In these models, there are also singlet scalars S within the same supermultiplet as S. 
These very weakly coupled scalars have masses of order mw (or, in the case of the invisible 
axion, essentially 0). However, decays like NiSS depend on couplings that are effectively 
doubly suppressed by e. Other decays involving only singlet scalars always have competition 
from ordinary unsuppressed MSSM decays, and so are not relevant for colliders. 

The rest of this paper is organized as follows. In section 2, I examine some specific models 
which realize the idea outlined above. Section 3 discusses the couplings and mixings of the 
singlino/ axino, and the relevant decays. Some representative numerical results for decays 
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to the singlino are presented in section 4. Section 5 contains some concluding remarks. An 
Appendix contains complete formulas for decay widths, including the effects of arbitrary 
phases. 



2 Singlino masses and couplings in models with an 
intermediate-scale solution to the fi problem 

Let us consider several models which realize the basic idea outlined in the introduction. The 
magnitude of / ~ (Xi j2 ) is (up to dimensionless couplings) the geometric mean of the Planck 
scale and the weak scale in order to agree with eq. ( |1 . 1| ) . One way that this could happen 
is if the soft supersymmetry-breaking (mass) 2 of X% is driven negative at an intermediate 
scale. More generally, X\ and X2 correspond to a nearly-flat direction in the potential, so 
that dimensionless supersymmetry breaking terms involving Xi and X 2 will always favor a 
non-trivial minimum at an intermediate scale. The key things we want to show are that 
these models generically contain one or more singlet fermions which have electroweak scale 
(or smaller) masses, and couplings that are of order v/f (with possible enhancements) to 
the MSSM Higgs fields. 

For example, suppose that the superpotential contains, in addition to eq. Q1.2Q , a term 

W = -h*-X x Xl, (2.1) 
6M P 1 2 ' K J 

as in H . This fixes the PQ charges of the superfields, ensuring the presence of an invisible 
axion provided that no other terms break the PQ symmetry. Here Ax is a dimensionless 
coupling which is assumed to be of order unity. The supersymmetry breaking Lagrangian 
must then include 

- £ S USYbrcaking = m l\ X l\ 2 + ™l\X 2 \ 2 - -2L(XiXf + C.C.) (2.2) 

blvl p 

where hx is a mass parameter of order the electroweak scale and has been taken to be real 
and positive without loss of generality. A nontrivial global minimum will exist provided 
e.g. that m 2 is negative at the scale of the VEV. However, it is important to note that this 
is not necessary. The presence of a holomorphic coupling hx always favors spontaneous 
symmetry breaking at an intermediate scale. So VEVs for X\ and X 2 can arise from a 
negative squared mass, and/or an hx which is sufficiently large [[H], 0. In any case, {Xi) 
and (X 2 ) are of order (m^Mp) 1 / 2 ~ 10 10 GeV, which is naturally within the invisible axion 
window. 



7 



Let us parameterize the VEVs of X\ and X2 by an overall magnitude / and an angle 
0, so that (Xi) = /cos0 and (X 2 ) = /sin0. Now expand the fields around their VEVs to 
obtain low-energy supermultiplet degrees of freedom S%, S 2 : 

f X x \ _ f coa<j>\ . f coaO ain9\(S 1 \ (0 „x 

\X 2 J ~ Uin^y 1 + cos6)\S 2 )' 1 } 



By requiring that the superpotential masses for the fermions S\ and S 2 derived from eq. ( |2.1| ) 
are diagonal, one can solve for the mixing angle 9 in terms of the VEV angle 0, with the 
result 9 = 0/2. [This typically does not diagonalize the axion and other light scalar masses, 
and depends particularly on the choice of eq. (|2.1|).1 The resulting masses and couplings for 
Si and 5*2 are then found to be 

m~ =M! sin( K C0S( />_iV v / cos 0/2 _ sin 0/2 \ 

5 1 2M P S 01C S0 j ' ^ A cos0 sin0 )' 1 j 

^xP •,/,,-, x w / sin 0/2 cos0/2\ 
m ? — 7TT7 sin (cos 0+ 1); = -? : — r~ • l^- 5 ) 

52 2M P v 7 S2 / y cos0 sm0 J y ' 

[Compare eq. (|1.3|).] The scale / and the angle could also be computed, in principle, in 
terms of the parameters in the soft supersymmetry-breaking Lagrangian eq. ( |2.2|) and the 
superpotential. The same parameters also determine the soft scalar mass of the saxion and 
other scalars with electroweak scale masses. However, these will not play any direct role in 
this paper, so I will not do this explicitly, and I will treat / and as free parameters. 

One interesting limit is that of small (i.e., (X 2 ) small compared to (Xi)), in which the 
mass eigenstate Si is the axino. Then one finds that, up to phases, 

m ^ = 4Mp- Sm0; ^ = 27 ; (2 - 6) 

^xP . , v 

rria = — — sm0: = -— — -. 2.7) 

& M P r ' 52 /sin y 1 

Note that both masses become small in this limit. The coupling eg- must grow like 1/ sin0 
in this parameterization in order for /i to not become much less than v. (LEP has not 
discovered a higgsino.) Another interesting limit is = 7r/4 (VEVs of equal magnitude), 
resulting in 

m * = o - io %T ; e?i = °' 77v/f; (2,8) 

m~ S2 = 0.60^; eg, = 1.85v/f, (2.9) 



S 



again, up to phases. Both Si and S% have masses that are roughly of order the electroweak 
scale. In general they each contain an admixture of the axino. It is interesting to note that 
\ m s 1 \/{^xf 2 /Mp) < 0.1 over the whole range 6 < tt/A. So it is not unlikely that one or both 
of S\ and 5*2 is lighter than all MSSM sparticles. 

Other models can be obtained by choosing superpotentials 

Ax ~q A, 



as in ref. 22 



or 



W 



w 



6M P 2 2M P 2 ' 



Xx -XiXl+ ^ ^ 



X 1 H u H d 



(2.10) 



(2.11) 



6M P * 2M P 

In both of these cases, the diagonalized singlino masses are still given as in eqs. ( |2.4|) and 
(p.5|) , since they only depend on the Ax term in the superpotential. However, the couplings 
are modified to, respectively: and 



-Si 



for eq. (|2.10|) , and 



Si 



2v sin cf)/2 
f smcj) 

2v cos 0/2 

/ COS (j) 



■s 2 



2v cos (j>/2 
f sin cf) 

2v sin 0/2 
/cos0 



(2.12) 



(2.13) 



for eq. (|2~TT| ). 

Another similar model is obtained by assuming a different form of the Ax term used to 
stabilize the potential at large field strengths: 



W 



A; 



-x\xl + 



A, 



-X 1 H u H d . 



4Mp" L " z ' 2Mp" L ~"""' ^ 2 ' U " > 
Again using the mixing parameterization eq. ( |2.3j ), one finds in this class of models that now 
6 = (1/2) tan -1 (2 tan 20) in order that Si and S2 are mass eigenstates. In terms of the VEV 
angle 0, we find: 

Ax/ 2 



m 



Si 



m 



s 2 



s 2 



4Mp 

V2v 
/cose 

Ax/ 2 
4M P 

V2v 
/cose 



1 - cos 20 1/1 + 4 tan 2 20 



1 + 



1 



1 + 4 tan 2 20 



1/2 



1 + COS20V1 +4 tan 2 20 



1 - 



1 



1 + 4 tan 2 20 



1/2 



(2.15) 
(2.16) 
(2.17) 
(2.18) 
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up to phases. In the limit — > 0, one finds 

m ^ = Tif sin20; e ^ = T (2,19) 

m-s = - - - ; = —smd>, (2.20) 

s 2 2M P 2 f 

which features a parametric suppression in the mass of Si. In the case of equal VEVs 
<f) = it /A, one obtains instead 

(2.21) 
(2.22) 



Xxf 2 


2w 






3Ax/ 2 









Finally, the limit — > n/2 (i.e., (Xi) small compared to (X 2 )) yields a very light S 2 and a 
relative enhancement of the Si coupling: 



Xxf 2 2v 
m ^ = 2Mp 5 ^ = /W ^ 



m ~ s - = ^Mp~ cos2 0; e " 2 = 7' (2 ' 24) 

There are clearly many possible more complicated variations on these models; for exam- 
ple, schemes with more than two fields X, participating in the PQ-breaking and /^-generating 
dynamics (see, for example, ref. |2c|). It is also possible to have a scheme in which there is 
only one field X, which obtains a VEV at an intermediate scale below where the soft mass 
term m 2 x runs negative. This corresponds to the <fi — > 0, Ax — > limit of eq. ( 2.14j) with 



S 2 removed, so there is just a light axino with e = 2v/f. The essential features of all these 
models are that they contain one or more singlino fields, with couplings naively of order 
v/f but which can be significantly enhanced, and which can easily be lighter than all of the 
MSSM odd--R-parity sparticles. 



3 Mixing of the singlino with MSSM neutralinos 

As shown in the previous section, one or both singlino mass eigenstates Si or S 2 can be 
lighter than all MSSM sparticles. In this section, I will consider the relevant mixings and 
couplings of such a singlino to the MSSM states, and the ensuing decay partial widths. I 
will use S to refer generically to either Si or S 2 . 
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The properties of the singlino are determined by the superpotential eq. ( |1.3[) . At tree 
level, there are singlino-higgsino-Higgs boson couplings. Other couplings, including singlino- 
fermion-sfermion and singlino-higgsino-Z boson, arise due to singlino-gaugino and singlino- 
higgsino mixing. In order to discover the couplings of the S to the physical MSSM states, 
one must diagonalize the the 5x5 neutralino mass matrix. In the (S, B, W°, H^, H®) basis, 
it is given by: 










Mi 


-c l3 s w m z 
s l3 s w m z 






M 2 
Cp c w m z 
-S/3 c w m z 



-c l3 s w m z 
op c w m z 


-fi 



-e/icp \ 
Sps w m z 
-sp c w m z 
-/i 

/ 



(3.1; 



where S/3,cp stand for sin/?, cos/?, and sw,cw for sm9\y, cos^iy. In the following, I shall 
take nig to be real and positive without loss of generality. This allows e to have an arbitrary 
phase. Now, the off-diagonal terms proportional to e can be treated as a perturbation. 
Therefore, our procedure is to first diagonalize M«, the lower right 4x4 mass sub-matrix. 
This is accomplished with a unitary matrix Z^ (i, j = 1, . . . , 4) according to: 



Z* k M£>Z. 



(4), 
kl ^jl 



(3.2) 



Here the masses are real and positive; this can always be done, regardless of the relative 
complex phases of /i, Mi and M 2 . To the lowest order in a perturbative expansion in e, the 
singlino S is a mass eigenstate, and the ordinary MSSM neutralinos have the same masses 
that they would have had if S were absent. I will choose an ordering scheme such that 
S = N , with = < < < < m^. 

The full 5x5 neutralino-singlino mass matrix can then be diagonalized according to 



(3.3) 



where now i, j = 0, 1, . . . , 4. To lowest order in a perturbation in e, one finds that iVoo = 1 
and Nj. 



Zij for i, j 



1, 2, 3, 4. So one can write 



1 

N l0 



N 0j 



The neutralino-singlino mixing elements can be determined in terms of the 4x4 Z^, 
mass eigenvalues, and the parameter e: 

4 



(3.4) 
the 



N< 



o.i 



J2 Z kj e*fi*m^ k (spZ k3 + c p Z kA ) + efim^spZ^ + c p Z* kA ) / (m 2 ~ u - m 



k=l 



2 



N, 



iO 



[e//m^. (spZ* 3 + cpZ* A ) + e*/i*m^(s /3 Z i3 + cpZ u ) / (m 



4, 



(3.5) 
(3.6) 
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(Again, nothing has been assumed here about the complex phases of the parameters M±, 
M 2 , [i and e.) This expansion is an excellent approximation because |e| is very small. In 
particular, one can check by numerical diagonalization of eq. (|3.1|) that the perturbative 
result is very accurate unless the denominator m~ — m~ in eqs. ( j3i>D and fl3.6|) is tuned to 

to an extreme accuracy (comparable to |e^mgj), in which case the decay widths studied 
below are completely negligible anyway. 

The relevant neutralino and singlino couplings for low energy phenomenology are all 
contained in the five eigenmasses and (for which corrections proportional to e are 
negligible) and the neutralino mixing quantities AV,- and Nqj with = 1,...,4). The 
decay rates with a singlino in the final state can now be worked out, as in the NMSSM Jl2| . 

1 present formulas for the relevant widths in an Appendix, taking care to allow for possible 
arbitrary phases in the neutralino mixing matrix and in e. In general, the phase of e could 
be anything, since it is not constrained by low-energy experiments on CP violation and is 
not necessarily correlated with the phase of \x. Because the results might be useful for other 
problems, I will provide the general results for a decay involving Ni to A^-; the case needed 
in this paper is obtained by simply taking i — 1 and j = 0. Numerical illustrations of these 
results will be given in Section 4. 

The necessary coupling for the decay 

iVx -> h°S (3.7) 

arises directly from the superpotential eq. ( |1.3| ) through the higgsino content of N%. It also 
obtains contributions from the singlino-higgsino mixings N 03 and A r 4 combined with the 
gaugino content of N\, and from the singlino-gaugino mixing elements Aoi and A/02 combined 
with the higgsino content of Ni. The total coupling is given explicitly by eq. (|A.2 ). The 



Higgs will then decay according to h° — > bb, WW*, t + t~, cc, or gg, providing a visible 
product displaced from the original interaction point producing the event. 
The decay 

Ni -> ZS (3.8) 

relies both on singlino-higgsino mixings A03 and A r o4 and on the higgsino content of Ni, 
and is therefore somewhat more suppressed in models with a bino-like NLSP. The relevant 
coupling is given explicitly by eq. (|A.6|) . The Z boson then decays with Standard Model 
branching fractions to quark-antiquark and lepton-antilepton pairs. 

When neither two-body decay is open, the neutralino NLSP will decay through off-shell 
sleptons, squarks, the Z, and Higgs bosons. The full expressions for these decays are given 
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in the Appendix. The most important contribution for a bino-like NLSP is typically through 
sleptons, and therefore relies on the singlino-bino and singlino-wino mixing elements N i and 
iVo2- Fortunately, these are not greatly suppressed in many models unless \fj,\ is very large. 

There remains the possibility of a two-body decay Ni — > jS, which arises at the one- 
loop level. However, using the results of [^4] (with appropriate changes for the couplings 
to correspond to the model under present consideration), I have verified that these decays 
are always suppressed compared to the two- and three-body decays considered here, with 
widths that cannot be much larger than a few times (1000 km) -1 in the examples in the 
next section with e = 10~ 8 . Since these decays are not competitive here, I will not present 
results for them. 

A separate possibility is that the NLSP is a stau, or that all three lighter, mostly-right- 
handed, slepton mass eigenstates (ri,^, e^) have no open decays except to the singlino. In 
that case, one can hope to observe f\ — > rS (and perhaps p,R — > fiS and en — ► eS). The 
decaying slepton will appear in the detector as a muon-like charged particle track, or as a 
track with an anomalously high ionization rate. The rare decay to the singlino will yield a 
large-angle kink in the track leading either to a tau jet or an electron or muon. Since the 
decaying particle is heavy, there will be a significant angle at the kink. The decay widths 
are suppressed only by the singlino-gaugino mixing, so they can occur within the detector 
often enough to measure, even if e is significantly less than 10~ 8 . 

4 Representative results for decays to the singlino 

In this section, I will consider some illustrative numerical results for decays to the singlino, 
first for neutralino NLSP models and then for stau or slepton NLSP models. I will take 
e = 10~ 8 , with the understanding that the results have to be scaled according to T oc e 2 . 

4.1 Neutralino decays 

In order to study the decay partial widths of a neutralino, I will employ the concept of 
"model lines", in which one supersymmetry-breaking parameter is allowed to vary, setting 
the overall scale for all sparticle masses. First, consider a typical model scenario with a 
bino-like NLSP and the LSP singlino mass fixed at m^r = 50 GeV. The bino mass parameter 
Mi is varied, with the wino mass parameter M 2 and the \i term then determined according 
to M2 = 2.0Mi; fi = 3.0Mi. The 5x5 neutralino mass matrix is then fully determined by 
also choosing fixed values of tan/? = 3.0 and e = 10~ 8 . The right-handed slepton masses 
m~ = m~ = m~ are constrained to be the greater of 1.2m^ and 110 GeV. This assures 
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10° 




N 1 mass [GeV] 

Figure 1: Visible decay widths of bino-like neutralinos N\ into final states involving the 
singlino, as a function of varying for fixed = 50 GeV and e = 10~ 8 . The MSSM 
model parameters are described in the text. The solid line is the total visible decay width. 
The long dashed line is the partial width into £ + £~S where £ = e or //. The dot-dashed 
line is the width into t + t~S. The short-dashed line includes jjS where j is any u,d,s,c 
(anti)-quark jet or gluon jet, and the dotted line is for bbS. The thin solid line is the width 
for jyW'iS through an on-shell Higgs boson. 

that a slepton cannot be the NLSP and should not be found at LEP. Mixing in the stau 
sector is neglected. The left-handed slepton masses are determined by m~ = m~ +0.5M|. I 
will assume that squarks are not light enough to give a significant contribution to the decay. 
Finally, the lightest Higgs boson mass is assumed to be m^o = 120 GeV, safely out of the 
reach of LEP, and to obey the decoupling limit a = /3— tt/2. The results for the partial decay 
widths to visible states (excluding neutrinos), as found from the equations in the Appendix, 
are shown in Figure [j] SIS db function of mr r . 

For rrifi < 120 GeV in this model line, the decays are dominated by the contributions 
of the virtual right-handed sleptons. The total inverse decay lengths are of order tens of 
kilometers, and are nearly democratic between e + e~, [i + fi~ and r + r~ final states. The 
"knee" near = 92 GeV is merely an artifact of the constraint m~ fl > 110 GeV; for smaller 
neutralino masses, the virtual slepton is necessarily more off-shell due to the constraint that 
it has not been discovered at LEP. 
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For > 120 GeV, the contributions from the virtual Z boson [the terms Wz, Wzu 
and Wzu m ec L- (|A.12]) 1 begin to be important. This increases the decay partial widths, with 
contributions to // that are roughly proportional to the Z branching fraction, so that dijet 
final states dominate. For > 142 GeV, the virtual Z boson is on-shell, and the decay 
becomes two-body Ni — ► Z°S. (The three-body formula is used in the vicinity of threshold, 
however, in order to correctly include interference effects with the virtual slepton diagrams 
in that regime.) This leads to a total visible decay width greater than (1000 meters) -1 . For 
negative /i, the decay widths tend to be somewhat smaller. 

Finally, for > + m^o = 170 GeV in this model line, the decay N\ — > h°S opens 
up and completely dominates. Since there is a direct higgsino-singlino-Higgs boson coupling, 
this is much larger than the two-body decay to ZS, even though the Z boson is lighter. The 
results are shown assuming Standard Model branching fractions for h° into final states bb, 
WW*, r + r - , and (lumped together into the "jj" category) cc and gg. Here the partial 
decay width of N\ to the bbS final state is found to be of order (10 meters) -1 . Of course, if 
the h° mass is smaller, this mode will open up and dominate for smaller values of . 

In the era after supersymmetry is discovered, the situation will be rather different; we will 
presumably know the MSSM sparticle mass spectrum, but the singlino mass and coupling 
will be completely unknown. So, a more useful summary of the situation we could face 
might be something like that shown in Figure |2|. This is a particular point along the same 
model line just discussed, with = 150 GeV, and with the horizontal axis representing 
the possible values of the singlino mass. 

A somewhat different scenario ensues if the NLSP is a higgsino-like neutralino. To illus- 
trate this, I choose a pair of model lines with fi = ±0.8M 1; and all other parameter relation- 
ships as described above for Figure [1]. The results are shown in Figure |3], but now only for 
the total visible decay width. Since Ni has a smaller gaugino content, the decays through 
virtual sleptons are highly suppressed. Conversely, the large higgsino component of N\ en- 
hances the probability of decay through a virtual Z boson. So, for < nig + m^ = 170 
GeV in this model line, the // decays will obey Z boson branching fractions. However, there 
turns out to be an accidental suppression of the matrix element for \x < in the convention 
specified in eq. ( [3.1|) (which is the same as in refs. [[I], [KJ), particularly when the Z boson 
is off-shell. For > 170 GeV, the decay Ni —>■ h°S length is of order several meters if Ni 
is mostly higgsino. 

In the above analyses, I have assumed that H° and A are very heavy, and that three- 
body amplitudes involving them are negligible. Although this is appropriate throughout 
most of parameter space, it is possible that for large tan j3, the couplings of H° and A could 
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LSP mass [GeV] 

Figure 2: As in Figure 1, but for fixed uit, = 150 GeV, and varying S (LSP) mass. 




50 75 100 125 150 175 200 
N 1 mass [GeV] 



Figure 3: Total visible decay widths of higgsino-like neutralinos N% into final states ffS, 
as a function of varying for fixed = 50 GeV. The model parameters satisfy the 
constraints M 2 = 2.0M 1 , tan/3 = 3.0, and fi = +O.8M1 (solid line) or fi = -0.8M l (dashed 
line), and other constraints described in the text. 
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be large enough to make an appreciable contribution to T(N\ — > bbS), even if H° and A 
are far off-shell. 

Finally, I note that the decay widths discussed here depend on the phase of the parameter 
e. This phase is not constrained by low-energy CP violating observables, and so might be 
considered completely arbitrary. I have checked, using the formulas in the Appendix, that 
varying Arg(e), while keeping all other parameters fixed, can change the Aq decay widths 
by an order of magnitude or so. In models with a bino-like NLSP as in Figure 1, the largest 
decay widths tend to occur for real e. 



4.2 Slepton decays 

It is also possible that the NLSP is a slepton fi or, effectively, all three mainly right-handed 
sleptons eR, JIr and f\. The latter scenario is realized if tan/5 is not too large, so that the 
three sleptons are mass-degenerate to within less than m T . These possibilities are familiar 
^-[fTi] in gauge-mediated supersymmetry breaking models |^2| but could also be realized in 
supergravity-mediated models if there is not a large universal contribution to scalar masses, 
or if .D-term contributions (proportional to some exotic £7(1) quantum number) are large. 

If T\ is the NLSP, then the two-body decays T\ — > rS are suppressed only by the bino- 
singlino mixing. In Figure f|, I show the results for this decay width for a typical model 
line with varying Mi and fixed LSP mass m^. In order to ensure that a stau is the NLSP, 
the constraint m~ = 0.9m^ (thick lines) and 0.7m^ (thin lines) are imposed. Other 
relevant model line parameters are M2 = 2.0Mi, /1 = 3.0Mi, tan/? = 3.0. To a good 
approximation, the decay width depends only on the absolute value of s~ in the stau mixing 
parameterization of eq. ( |A.24|) , so results are shown for s~ = 0, 0.25, and 0.5. The solid line 
is also approximately true for eR — > eS and JIr — > fiS as a function of m~ R and rn~ , by 
taking s~ = 0, c~ = 1. 

Because there is relatively little suppression in this case, the inverse decay widths are of 
order tens of meters. This is discernible at a collider which can produce several hundred 
super symmetric events. Each such event would contain a pair of quasi-stable stau or slepton 
highly ionizing tracks, which can have an anomalously high dE/dx to distinguish them from 
muons. In a small fraction of events, one of the stau or slepton tracks will have a kink leading 
to a lepton or tau jet, corresponding to the decay. The resulting tau or lepton would have 
a significant angle with respect to the original highly ionizing track, yielding a potentially 
spectacular and nearly background-free signal. 
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Figure 4: Decay widths for f\ — > tS as a function of varying m~ ■ . The model parameters 
satisfy mg- = 80 GeV, and the constraints M 2 = 2.0Mi, = 3.0Mi, and tan/5 = 3.0. Solid 
lines are for the unmixed case s~ = 0, the dashed lines for s~ = 0.25 and dash-dotted lines 
for s~ = 0.5. The thicker (thinner) lines are for m~ /mj^ = 0.9 (0.7). 

5 Conclusions 

The presence of the \i term in the MSSM and the solution to the strong CP problem may have 
a common explanation at an intermediate scale. Direct detection of the resulting axion is 
quite problematic. In this paper, I have argued that these models may nevertheless give rise 
to observable signals at colliders, through delayed decays to singlino fermions that include 
the axino as a mixture. These events will be a rare (perhaps very rare) occurrence within a 
large sample of supersymmetric events at the Large Hadron Collider or a future e + e~ linear 
collider. There is also a possibility that the lightest MSSM sparticle could have slow decays 
into more than one singlino. Note that the LHC cross sections can be very large precisely 
when the NLSP decay widths are small. 

The numerical estimate in section 4 of this paper have used e = 10~ 8 for the singlino- 
higgsino-Higgs coupling parameter. Of course, the actual value could be significantly smaller. 
On the other hand, I showed that in some models the couplings are parametrically enhanced, 
and the mass of one or more singlinos is reduced, if one of the VEVs giving rise to the \x term 
is relatively small. Furthermore, the magnitude of e can be significantly larger if the high 
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scale Mp is replaced by a somewhat lower scale that governs non-renormalizable operators, 
for example a string scale or a compactification scale that is not far above the apparent 
gauge coupling unification scale. 

Future planning and analysis of collider physics experiments should take into account the 
possibility that the apparent LSP is actually unstable. Besides the models I have discussed 
here, there are at least two other plausible variations on the MSSM which can lead to delayed 
rare decays of what might appear, at first, to be the stable LSP. 

First, gauge-mediated supersymmetry-breaking (GMSB) models |32| with a supersymmetry- 
breaking scale y~F that is not too large will give rise to decays that could have macroscopic 
proper lengths ||25||-||31j|. It is interesting to compare the reason for this to that in the models 
discussed in the present paper. In GMSB models, an estimate for a decay width of the NLSP 
to the goldstino/gravitino G is V ~ m^/167r(v / ^) 4 , while in the decays to a singlino/axino 
LSP S, the estimate is T ~ / 167T f 2 . So NLSP decays are suppressed by the 4th power 
of the supersymmetry-breaking scale in GMSB, but only by the square of the PQ scale in 
light axino/singlino models. GMSB models can, in fact, give rise to signals which might 
be very difficult to distinguish from those discussed here. For example, if the NLSP is a 
neutralino with a significant higgsino content, it can have |53L |2B|, |34J decays Ni — > h°G 
and Ni — > Z°G that look like the decays discussed in this paper. Or, if a stau is the NLSP, 
it can appear quasi-stable with rare decays fi — > tG occurring within the detector. Second, 
one can have weak i?-parity violating couplings in the MSSM which could also give 
decays like Ni — > t + t'~v or N\ — > qq'v. These signatures could mimic those discussed in the 
present paper. 

If these signals appear, it will be interesting to try to establish the correct explanation 
from among the competing hypotheses. The prize for doing so will be that we will gain an 
understanding of physics at scales far above those probed by direct sparticle production at 
colliders. 
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PHY-9970691. 



Appendix: Complete formulas for neutralino decay widths 
including the effects of arbitrary phases 

In this appendix, I give formulas for the two- and three-body decays of a neutralino to 
another neutralino and a Higgs boson, Z boson, or fermion-antifermion pair. The model is 
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the extension of the MSSM with one singlet superfield as specified by eq. ( |1.3|) . This also 
includes the MSSM and NMSSM as special cases. The results needed in Section 4 of this 
paper are obtained by taking the neutralino mass eigenstate indices to be i — 1 and j = 
in the following. 

First, let us consider the two-body decay of a neutralino to another neutralino and the 
lightest CP-even neutral Higgs boson. The relevant decay width is equal to: 

T{%^h%) = ^| V /A(l,r],r2 )(|Gg , | 2 (l +rj - r 2 ) + 2Re[(Gj°) 2 ]r J ) (A.l) 

where Tj = ; r h o = m h o/m^ ; A(a, b, c) = a 2 + b 2 + c 2 — lab — lac — 2bc; and the 

neutralino-neutralino-Higgs coupling is given by 



G<S = ^e-^X*^^ (A.2) 



Here c a and s a denote cos a and sin a with a the Higgs mixing angle in the notation of [10 



The results for the decays to the heavier neutral CP-even (H°) and CP-odd {A ) Higgs 
bosons can be obtained by instead using the couplings: 

Gf- = ^*-s'iV*J(- c ^ (A.3) 
for % -> H°Nj, and 

Gf- = ~(^*-s'Ar*)( s ^^ (A.4) 

for Ni — > A°Nj, and substituting m^o — > m^o or m^o in the obvious way. However, in the 
numerical analyses of Section 4, H° and A are assumed to be heavy and decoupled, so these 
decays are neglected. 

Similarly, the two-body decay of a neutralino to another neutralino and a Z boson has a 
width given by 

r(iV, Z%) = ^ v /A(l,r 2 ,r 2 )(|G|| 2 [l + r) - 2r| + (1 - r 2 ) 2 /r|] + 6Re[(Gg) 2 h) (A 
where = m^/m^ , and the neutralino- neutralino- Z coupling is given by 

G% = 2^(-WV*s + ^iV; 4 ). (A.6) 

Finally we consider three-body decays iVj — > ffNj, where / is any Standard Model quark 
or lepton. Results for these decays have appeared in refs. [O, |36], 571, but here I include 
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the effects of Higgs boson exchanges and arbitrary phases of the couplings . The differential 
partial widths can be expressed in terms of the dimensionless mass ratios r/ = mf/m^ , 
Tj = mfijmfr, r z = m z /m^ , rr z = T z /m^ , and = m^/m^ for each of = h°, A , and 
H° . I use dimensionless kinematic variables 



N, 



U 



1 + r] + 2r 2 f - s - i 



with limits of integration 



1 



^ mm, max 



= 4rJ; 



2 [1 + rj - s + 2rj T {(1 - 4rj/s) A(l, s, rj)} 

^max (1 fj) ■ 



2 U 1/21 



The results for widths can be expressed as* 

driN, -> ffNj) = (£ w) dids , 

where n c = 1 (3) for leptons (quarks). The individual contributions to W are: 



W z = 



W t 



w u = 



M + 6 i 



(s — rV) 2 + r|rp. 



■{IG 



Z|2 



l-fi)(fi-rj) + (l-f)(f-r# 
+2Re[(Gg)>^} 



( r 2 _ t )( r 2 _t) 



tu 



2 n_fuf_ r 2) 

X + b]b]'*)(aT4 + !>■■■!>;')■ ' ' ' 

n,n'=l 

w t (i -> «) 

2Re £ 



W h o jH o 



4(s - r|) 



(a^ n a™ %T + a™ 6™ar&™ *) W - tu) 
+ (a]a]'a r l*a r ?'* + b]bfb™b?*)s rj 



(s-r 2 z ) 2 + r 2 z rf z „ , 
W zt (*->«) 



Re £ [(a,ay? + b f ^){G^\ - t)(t - rj) 

+G?*sr j }]/(rl-t) 



E 



=h°,H° 



4sRe[GjGj 



(A.7) 
(A.8) 
(A.9) 



(A.10) 
(A.H) 

(A.12) 



(A.13) 
(A.14) 

(A.15) 
(A.16) 



(A.17) 
(A.18) 



1 + rj - S)Re[Gg.Gg*] + 25r,-Ite[Gj.Gg]} (A.19) 



tin computing these results, I have neglected fcrmion masses arising from spinor algebra in matrix ele- 
ments, but not in the kinematic limits of integration or the couplings. 
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As\G A °\ 2 

W A o = pi ± J )2 {(1 + r) - 3)|Gffi| 2 + 2gr J Re[(G^) 2 ]} (A.20) 

<j>=h°,H°,A° n=l,2 l'<£ ^ /„ ' 

W*, = WV(*->«)- (A.22) 
Quantities appearing in the above results are as follows. First, 

a f = — — (T 3/ - qfSw)] b f = -q f gs 2 w /c w (A.23) 

Cw 

are the Z boson couplings to quarks and leptons with (T 3 f,qf) = (1/2,2/3) for up-type 
quarks, (—1/2,-1/3) for down-type quarks, and (—1/2,-1) for charged leptons. Left-right 
mixing and CP violation in the sfermion sector are parameterized by a unitary matrix, which 
I choose| to write as 



fR\ - ( c f s Af fj 



(A.24) 



where \cj\ 2 + |s^| 2 = 1, and mj^ < m j 2 , The resulting couplings for down-type sfermions 
(b, f) are: 

a} = s/^sjignfN^ + gXqf-T^N^-cjgN^mf/^mw^ (A.25) 
bf 1 = V2cjg'q f N il + sygN i3 m f /V2cpm w (A.26) 
a{ 2 = -^cjigTsfN^ + gXqf-Ts^N^-spN^mf/V^cpmw^ (A.27) 
b{ 2 = yfisjg'qfNa-cjgNamf/yficpmw (A.28) 

for i — 0, 1, . . . 4. For up- type fermions one must replace Nj£'/cp by /sp in eqs. ([A.25Q - 

( |A.28| ). (However, for the cases of interest in this paper, N± — *> tiS is surely not kinematically 

allowed, and decays Ni — > vvS are not interesting.) In the expressions for the contributions 

to the widths, I have used the abbreviations a™ = a{ n , etc. The various Higgs boson couplings 

to Standard Model fermions are given by, e.g.: 

h o _ gm b s a , gm t c a 

zmwcp 2m w S/3 
H o _ gm b c a H o _ gm t s a 

2m w cp 2m w sp 

G b = -i— ; G t = -i— , (A.31 

2m w 2mw 



■'"This parameterization has the feature that in the typical unmixed, CP-conserving case cj= b s j — 0, 
fl = fa an< i ]i = Sl, with no minus signs. 
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for bottom and top quarks. The Higgs couplings for taus are obtained by mj — > m r , and for 
other quarks and leptons by the obvious substitutions. 

The two-body decay width for a stau to a neutralino or singlino is given by: 

r(ri-TJVi) = ^ V / A(l,r i 2 ,r2){(|a[ 1 | 2 + \bf\ 2 )(l - r] - r 2 T ) - ^Re^^l} (A.32) 
where now = m^Jm- and r T = m T /m~ , and the couplings aj 1 , bj 1 are given already by 



eqs. (|A.25|) - (|A.28|) . In section 4, this formula is used with % — 0, corresponding to f\ — > rS. 



The results for [Ir —* and e# — > eS are obtained by taking c~ — > 1 and s~ — > 0. 
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